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ABSTRACT 

We report a search for z~6.96 Lya emitters (LAEs) using a Narrow-Band 
filter, centered at 9680A with the IMACS instrument on the Magellan telescope 
at Las Campanas Observatory. We obtain a sample of 6 Lya emitter candidates 
of luminosity ~ lO^^erg s~^ in a total area of 465 square arcmin corresponding to 
a comoving volume of ~ 72000M pc^. 

From this result, we derive a Lya luminosity function (LF) at z~6.96 and compare 
our s ample with the only z~6.96 Lya emitter spectroscopically confirmed to 



date ( live et al.l 120061 ) . We find no evolution between the z=5.7 and z~7 Lya 
luminosity functions, if a majority of our candidates are confirmed. Spectroscopic 
confirmation for this sample will enable more robust conclusions. 



1. Introduction 

Over the last decade, significant progress has been made in determining the processes 
of galaxy evolution using both ground- and space-based telescopes. Currently, the limits of 
the observable universe are at z ~ 6, which corresponds to ~ 90% the age of the universe. 
At z>6, we are approaching the NIR domain : the sky is brighter and it is more challenging 
to detect the faint high redshift sources. Beyond this boundary lie the first ultraviolet (UV)- 
emitting sources, which ionized the majority of the hydrogen in the universe. Their detection 
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will allow us to probe the era of reionization, after the "Dark Ages". Galaxies formed at 
high redshifts play a key role in understanding how and when the reionization of the universe 
took place. They also help constrain the physical mechanisms that drove the formation of 
the first stars and galaxies in the universe. 

Starbursting galaxies can emit a large fraction of their ultraviolet luminosity in the Lya 
line. Because Lya photons are resonantly scattered in neutral hydrogen, even a small amount 
of dust can quench this emission. Hence, selecting objects with strong Lyo; emission lines is 
expected to reveal a set of objects in the early phases of rapid star formation. These could 
either be young objects in their first burst of star formation or evolved galaxies undergoing 
a starburst due to a recent merger. Selecting galaxies with strong emission lines also allows 
us to probe the high-redshift Lya luminosity function (LF). 
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The presence of 

neutral hydrogen in the IGM can reduce the Lya flux of galaxies, it is therefore clear that the 
Lyo; LF is sensitive to the ionization fraction of the Universe. If we knew the intrinsic LF(z) of 
galaxies at each redshift, a deviation of the observed LF from this intrinsic distribution could 
be attributed to the attenuation by HI, and hence be used to infer the ionization fraction. 
In practice, the approach is to do proceed to a comparison of Lya LF at diffe rent redshifts, 



since the LFs of Lya emit ters (LAEs) don't evolve much between z=3-5.7 (ICassata et al. 
2010l : iMalhotra et al.l[201ll ). 



With ground-based telescopes, the detection of very distant objects requires observation 
of UV spectral signatures that have been redshifted into the visible spectrum. The longer the 
wavelength of the observed Lya line, the earlier the epoch at which we observe the galaxy, and 
the closer to the "Dark Ages" . Therefore, one way of searching for the most distant galaxies 
is to search for the redshifted Lya emission at the longest possible wavelength. However, this 
search is complicated by the presence of OH emission lines within the terrestrial atmosphere, 
at an altitude of ~ 80A;m. This strong line emission limits the sensitivity of ground-based 
telescopes at near-infrared wavelengths. Fortunately, there are spectral intervals with lower 
OH-background that allow for a fainter detection limit from the ground. 

We use a custom-built filter, A^i?9680, centered at A = 9680 A and with a width of 
QOAto use one of the low-sky windows. Narrow-band imaging is the most successful method 
to detect strong Lya emission lines of galaxies, since it relies on a specific redshift interval 
as well as a selected low-sky background window. Adapting the spectral width of this filter 
allows for maximum detection of light from the celestial objects at that spectral line, while 
minimizing the adverse influences of sky emission. 
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The first z > 6 LAEs was d etected with the narrow-band (NB) technique at the 10m 
Keckll telescope (IHu et al.ll2002l ). This galaxy was spectroscopically confirmed to be at z = 
6.56. Over 1,000 LAEs have been photometrically selected and spectroscopically identified 
in this way. Extensive observations have been done at redshifts 5.7 and 6.5, two spectral 
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confirmation yet. 

We present here a new NB imaging survey with the IMACS /Magellan telescope - targeting 
z = 6.96 LAEs. This paper first presents the data (Section 2.1) and the data reduction 
procedure (Section 2.2). We then describe the method of selection and contamination of 
low-redshift interlopers for high redshift LAEs in Section 3. We present the final sample of 
z = 6.96 LAEs and Lya luminosity function at this redshift in Section 4. 



Throughout this study, we adopt the following cosmological parameters : Hq 
= 0.3, Qa = 0.7 flSpergel et al.ll2007f ). All magnitudes are AB magnitudes. 
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2. Observations and Data Reduction 
2.1. Observations 

The data were taken with the IMACS instrument (the Inamori-Magellan Area Camera 
& Spectrograph), installed at the 6.5m Magellan Baade telescope at Las Campanas Obser- 
vatory. This instrument offers two cameras: f/2 and f/4, with different imaging scales and 
dispersions. We observed with the f/2 camera which delivers an image of 27.4' diameter field 
at a scale of 0.2 arcsec per pixel. To ensure background limited performance, each exposure 
lasted 15 min. After each exposure, the telescope was offset and a new exposure was taken. 
The offset values, for each exposure, was randomly chosen. 

We targeted one COSMOS field (RA=10:00:29 Dec=02: 12:21) covered by the Canada France 
Hawaii Telescope Legacy Survey (CFHTLS) and the WIRCam Deep Survey (WIRDS : Pis 
Willott & Kneib)|f|. The total area of the survey is 572 square arcminutes. We observed 
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using the NB filter centered at 9680A (A^59680) on the 17th, 18th and 19th March 2009. 
During these observations, the conditions were good : we obtained 19, 23, and 19 15-minute 
exposures during these nights, with a seeing varying between 0.45" and 0.8". 

We observed the field again on the 21st, 22nd and 23rd March 2010. Between the 2009 
and 2010 observing season, the filter was replaced with a new one built to have identical 
bandpass specifications. There is no major difference between the two filters' realized band- 
pass, and we treat data from the two identically. The night conditions were not as good 
as those in 2009 : the seeing varied from 0.7" to 1.3" and we lost on average 2 hours per 
night due to instrument problems. We obtained 16, 19, and 18 15-minute exposures in three 
nights. The total exposure time for each epoch of data is given in Table [H 

2.2. NB Data Reduction 

2.2.1. Data reduction 

To process this data, we use the package MSCRED/IRAF. At the beginning of each 
night, we take 10 bias frames and 10 dome flats, we process the bias frames by overscan 
subtraction, trimming, and stacking to produce a master bias frame for each night. The 
science frames are bias subtracted and flat fleld corrected in the standard way. But the main 
difficulty was the fringing correction. 

We see fringes produced by the interference of light reflected between parallel surfaces 
in an instrument. All object frames do not always share the same fringe pattern because 
flexure and variations in illumination geometry can change its amplitude or period even on 
short timescales. In order to correct each individual image for fringes, we produce a median 
image per night. An illumination frame is also generated for each science frame using a 
two-stage median smoothing algorithm, applied on each frame with a flrst stage flltering on 
a 16 pixel scale, the second on a 24 pixel scale. The net effect is comparable to a 384 pixel 
median fllter corresponding to 76.8 arcsec. 

We then produce a median of the illumination frame for each night and correct the 
science frames for the illumination pattern. We subtract this illumination frame from the 
median of the science frame to obtain a fringe pattern for each night. It is then essential 
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to find the optimal multiplicative scaling factor to use in removing the fringe pattern from 
each science frame, as we have long exposure frames and the fringe pattern changes with 
time. Once the frames are corrected for the fringes, we perform a sky subtraction by sub- 
tracting a normalized median image of all science frames from each frame. Finally, once 
the sky-subtracted frames are reconstructed as single extension images, we stack them using 
the mscstack /IRAF task. We made stacks per night, per epoch and a combined epoch stack. 

Astrometric calibration. We perform the astrometric calibration on the individual im- 
ages before the final stacking. We set an initial WCS information in the header of the frames 
based on the COSMOS catalog. We need then to adjust the WCS parameters to obtain 
a more precise alignment. For this purpose we use mscwcs to apply a first offset on the 
(RA,Dec) coordinates and we refine the calibration by checking the alignment with mscc- 
match task. We finally obtain an astrometry calibration for each individual images with a 
precision of rms~ ±0.1 arcsec in both directions. 



Photometric calibration. The photometric calibration of the CFHTLS data is based on 
the SP SS data for stars with 17<i'<21 and the Megacam-SDSS color transformation equa- 
tions of iRegnault et al.l (120091 ). The precision obtained in u* , g', r', i', z' is between 0.03 and 
0.02 mag. As the A^E9680 filter is included in the z'-band filter, we used the MAG AUTO 
magnitude from the z' band SExtractor catalog to calibrate our A^i?9680 filter. We select 
1048 non-saturated stars in the magnitude range 16< z' <20 to perform the calibration. 
Considering the photometric error on the broad-band (BB) calibration, we obtain a photo- 
metric calibration precise to 0.1 magnitude in NB9680. 



Completeness. We estimate the limiting magnitude for the different bands by adding 200 
artificial star-like objects per bin of 0.1 magnitude in blank regions of the different stacked 
images. We then run SExtractor on the image with the same parameters as previously used 
for object detection. We repeat this procedure 40 times. The average count on 40 times of 
the number of artificial stars retrieved in each magnitude bin provides a direct measure of 
the completeness limit. We report the 50% completeness limit in Tabled! 
We use this result to determine the Luminosity Function presented in Section |H 
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2.3. Broad Band Data 

The CFHT-LS provides extremely deep optical imaging data for our observed field. For 
the purpose of this study, we made use of the T0006 release. The CFHT-LS data products 
are available from the CADC archive to CFHT users and take form of image stacks in the 
u*, g', r', i', z' filters and of ancillary data such as weight maps, catalogs etc. The spectral 
curves of the filter u*, g', r', i', z' are similar to the ones of SDSS filters. In addition to the 
optical data, we also have the WIRDS survey, providing J, H, i^s. These optical data have 
been calibrat ed photometrically usin g the SDSS photometry and the NIR data using 2MASS 



photometry (IMcCracken et al.ll2010l ). Considering internal and external photometric error 
sources, the uncertainty on the optical and the NIR data photometry is ~0.05 mag and 
~0.02 mag, respectively. 

As the broad-band data and the narrow-band data do not have the same pixel scale, we 
resample the broad-band data using the software SWARP to obtain optical and NIR images 
with 0.2 arcsec/pixel. The alignment in pixels is then verified using IRAF geomap/geotran 
tasks. Our complete set of data is therefore scaled at 0.2 arcsec/pixel and covers an area 
of 27.4' diameter. A summary of the observational data used in this paper is provided in 
Table [H Figured] shows the transmission curves of the filters corresponding to the multi-band 
data used in this study. 



3. Sample 
3.1. Catalog generation and Selection 



We generate the catalogs using the software SExtractor (iBertin fc Arnoutslll996l ). We 
use the dual-image mode : the first image, for the detection, is settled as the combined 
A^i?9680 image, the second image, for measurement, corresponds to the resampled images 
from the optical and NIR bands. We choose to detect objects in 7 pixels above a threshold of 
1.5a, corresponding to a A^i?9680 ~ 25.8. The aperture used for the photometry is 1 arcsec 
diameter. 

As the IMACS instrument is composed of eight chips, we see an increase of the noise 
in the inter-chip regions. The sky noise is higher by a factor 2.5 in the interchip regions. 
We choose therefore to eliminate these regions from the area of the survey. From the total 
area covered by our survey, 572 square arcminutes, we obtain an effective area to search for 
z~6.96 LAE candidates of 465 square arcminutes. 
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Table 1: Observational data. 



Instrument 


Band 


Integration 


Limiting 






time (hrs) 


magnitude*^ 


MegaCam 


u* 


12.5 


26.6 


MegaCam 


9' 


20.25 


27.9 


MegaCam 




37.6 


27.6 


MegaCam 


i' 


54 


27.3 


MegaCam 


z' 


40.1 


26.4 


IMACS 


NB9680 1st epoch 


15.25 


25. 


IMACS 


iV59680 2nd epoch 


13.25 


25.*^ 


IMACS 


iV59680 combined 


28.3 


25.2^ 


WIRCam 


J 




24.8 


WIRCam 


H 




24.9 


WIRCam 


Ks 




25 



^ 5cr magnitude limits in apertures 1" in diameter for MegaCam and WIRCam. These 
limits correspond to a 50% completeness level. 

^ The filter used during 2010 observations has a better transmitted wavefront than the 
filter used in 2009. This explains the similar limiting magnitude between the two 
epochs although the observing conditions were different. 

By combining the two epochs of data, we should have expected a sensitivity increase 
of 0.37 mag. But the limit has increased by only 0.2 mag. This could be due to a 
systematic noise component like the fringes. 

I 




1000 



X(nm) 



Fig. 1. — Transmission curves of the filters corresponding to the data used in this paper. All 
transmissions include the response of the detector and are normalized to 100% at maximum. 
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Criterion^l : Since we have two epochs of data, we define A^i?9680 selection criteria 
on individual epochs and on the combined images. We select objects with a 3a detection 
in both of the individual epoch images and a 5a detection in the combined A^S9680 image. 
59% (=16150 objects) of the objects present in the initial catalog pass this criterion. This 
eliminates variable sources from the catalog. 



Criterion7^2 : Due to the nearly complete absorption of the flux shortward of Lya by 
the intergalactic Hydrogen, we should observe flux discontinuity at rest wavelength of 12 16 A, 
and observed wavelength of 9680A.We are therefore searching for objects which are not de- 
tectable in optical {u*, g', r', i') bands. A possible method is therefore to select objects with 
less than a 3cr detection in fllters blueward of the expected Lya emission : w*, r', i' . 
We use a image generated by combining g' , r' and i' CFHT data to obtain deep photom- 
etry of candidates in the combined optical bands. 3% of remaining candidates are accepted 
on this basis. 



Criterion^^S : Following the successful method of iRhoads fc Malhotral (120011 ) used 
for the search of z=5.7 LAEs in the LALA fleld, we require that more than 50% of the 
A^i?9680 flux comes from an emission line for object selection, which can be translated 
as magj^B9680 — t^clSz' < —O.Tbmag. This criterion is indicated in Figure |2l showing the 
z'—NB9680 color versus NB9680 magnitude diagram. 88% of the objects selected after the 
criterion#2, pass this color restriction. 



Criterion^4 : To avoid selecting extremely red objects (jCimatti et al.ll2002l ). we choose 
to set a different criterion : mag^Bgeso ~ i^o^dj < 0. 98% of the objects selected after the 
criterion^^S pass this color criterion. 



Criterion#5 : After applying carefully all the criteria presented above, we inspect each 
candidate visually. 33% of the objects inspected carefully are selected as serious candidates, 
the rest are rejected because they are near chip boundaries, or they are defects, etc. 

The equations summarizing this selection are presented in the Table [2j 
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Fig. 2. — Color-Magnitude Diagram z'-N B9680 vs NB9680 showing the candidates obtained 
with our criterion magNBdeso — t^agz' < —O.Tbmag (plain red triangle). The other points 
with magNBmm — rnagz' < —0.75mag are ruled out by other criteria (blue flux, variability 
etc.). 



3.2. Contaminants 



COSMOS photometric redsliift catalog gives us a first insiglit on tlie low redsliift emit- 
ters contaminating our high redshift candidates sample. 



Transient objects. Our strategy of observing the same field during two different years, 
and our requirement that the eligible candidates have to be detected in each epoch data 
within 3a level, allow us to rule out the contamination of our sample by transient objects 
such as supernovae, which would have appear in only one epoch of data. This strategy like- 
wise prevents contamination by slow-moving solar system objects. 



L-T dwarfs stars. Following the method described in iHibon et al.l (120 lOf ). we deter- 
mi ne the expected nu mber of L and T-dwarfs present in this survey. From the Figure 9 
of iTinney et al.l (120031 ) . representing the relation between the absolute magnitude and the 
spectral types of late-type dwarf galaxies, we find that we could detect L-dwarfs up to a 
distance of 871 to 3630pc and T-dwarfs up to a distance of 400 to 1260pc, depending on 
spectral type, from the coolest to the warmest. 

This field is located at high galactic latitude. Our sensitivity to L- and T-dwarfs is then 
extended beyond the scale height of the Galactic disk. H owever, the scale height applicable 
to L-, T-dwarfs is truncated at 350pc (IRyan et al.ll2005l ). We estimate therefore a sample 
volume of ~ 692pc^. Considering a volume density of L- and T-dwarfs of a few 10~^pc~^, we 
expect less than one L-, T-dwarf in our field. 

In addition, L-, T-dwarfs have A^i?9680-J>0 so they would have fail Criterion^4 of our se- 
lection. 



Lower redshift galaxies with strong emission lines should be actively star forming galax- 
ies with blue continuum emission. Looking for such a continuum allows us to identify these 



Table 2: Table of the selection criteria. 

Our criteria 
Criterion#l SNR{N B9680, ombined)> 5a 
SNR(A^S96802oo9)> Sa 
SNR(A^S96802oio)> 3(t 
Criterion#2 SNR(m*, g', r', i')< 3a 

Criterion7^3 nriagNBmm ~ i^CLg^' < —0.75mag 
Criterion 7^4 fnagNB9680 — ^o,gj < 
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lower-redshift line emitters, unless their equivalent qidth is very large. 

Foreground emitters. We estimate the minimum observed equivalent width a foreground line 

emitter would require to be selected with our criteria using the formula from lRhoads &: Malhotra 

fl200lh : 



EWmin ~ 



NB 



BB 



AX 



NB 



3aBB 



AX 



NB 



with /jVB and fsB the flux in A^i?9680 and g' band respectively, AX^b the width of the 
NB9680 filter, a^B and asB, the flux uncertainties in A^_B9680 and g' band respectively. 
We obtain therefore an EWmin ~ 1545A in observer frame. Foreground line emitters would 
then require an observed equivalent width EWmin ^ 1545A to contaminate our Lya selection 
sample. 

For guidance, this observed equivalent width corresponds to a rest-frame equivalent width 
of EWrest > 105lA for Ha emitters at z~0.47, of EWrest > 792A for [Olll] emitters at 
z~0.95 and of EWrest > 594A for [Oil] emitters at z~1.6. In the following studies, we used 
the observed equivalent width to estimate the number of emitters present in the survey and 
possibly contaminating our high redshift sample. 

1- Ha at z~0.47 

We flrst estimate the fraction of Ha emitters with EWobs > 1545A fr om Figure 2 of 
Straughn et al.l (120091 ). Fewer than 1.2% of Ha emitters at z~0.27 from IStraughn et al. 



(120091 ) sample would have such an equivalent width. We then evaluate the number of Ha 
emitters at z~0.47 present in our survey using the luminosity function from Figure 14 of 
Tresse et al.l (120021 ). We flnd that ~ 12 Ha emitters can be present in our A^i?9680 combined 
image. An upper limit of the number of Ha emi tters at z^O. 4 7 pas sing our criteria is then 
0.15. Considering the Ha luminosity function of iGeach et al.l (120101 ) . we flnd that ~ 52 Ha 
em itters can be presen t in our survey. An upper limit of the number of Ha emitters based 
on 



Geach et al.l (120101 ) is therefore 0.62. Thus Ha emitters are not serious contaminants. 



2- [Olll] at z-0.95 

We apply the same method for the Ha emitters to estimate the number of [Oi n j emi tters 
in our survey and contaminating our selection. From Figure 2 of IStraughn et al.l (120091 ) . we 
estimate the fraction of [ ill] emitters at EV Vrnin to be one out of 136. We obtain therefore 
an upper limit of 0.74%. iKakazu et al.l (120071 ) [O ill] emitter sample is covering a wide range 
of rest-frame equivalent width up to EWrest ~ lOOOA. As the rest-frame EW of the [O ill] 
emitters that could contaminate our high redshift sample is around 792.1 A we are therefore 
able to use their sample to estimate the possible number of [O ill] emitters passi n g thro ugh 
our selection criteria. From the luminosity function Figure 13 of iKakazu et al.l (120071 ). 48 
[O III] emitters could be present in our survey. Applying the upper limit to this number, we 
flnd that a maximum of 0.35 [Olll] emitters at z~0.95 could have high enough EW to be 
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selected. 



3- [On] at z~1.6 

We apply the same method to es timate the number of [O ll] emitters in our survey and 
contaminating our selection. From IStraughn et al.l (120091 ) , we obtain an upper lirn it of 3.3%. 
Using the luminosity function presented in Figure 5 of iRigopoulou et al.l fcOOSl ). 45 [Oil] 
emitters at z~1.6 can be detected in our survey. However a maximum of 1.5 of these objects 
can pass through our selection criteria. 



3.2.1. False Detections 

In order to estimate the number of false detections that could pass our selection criteria, 
we create an inverse A^i?9680 image by multiplying the combined A^S9680 image by -1. We 
then apply the same selection method and criteria and we did not find any candidates. 



3.2.2. Comparison with COSMOS redshifts 

We verify our sample of candidates by cross-correlating this catalog with the photomet- 
ric redshifts catalog from the COSMOqj field. This catalog covers a redshift range up to 
z~5.2 and a magnitude range from z' ~18 to z' ~25. None of the catalog's objects matches 
our candidate sample. Since most of the case against foreground emitters is discussed in the 
previous sections and a consistency verification performed, we can therefore conclude that 
it is very unlikely that our z~6.96 LAE candidate sample is contaminated by low-redshift 
interlopers. 

A more detailed study of the foreground emitters (Hce at ~0.47, [Olll] at ~ 0.95 and [Oil] 
at ~ 1.6) will be presented in a forthcoming paper. 



3.3. Final sample 

Our final sample contains 6 z~6.96 LAE candidates over the range of A^i?9680=24.1- 
24.4 and SNR (iVi?9680)=5.6-7.3, as described in Table Hi In order to derive a luminosity 



^This research has made use of the NASA/IPAC Infrared Science Archive, which is operated by the Jet 
Propulsion Laboratory, Cahfornia Institute of Technology, under contract with the National Aeronautics and 
Space Administration. 
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function independent of the photometry aperture, we compare the automatic aperture mag- 
nitude, the 1" aperture magnitude and the isophotal magnitude for unsaturated objects. 
We then correct the 1" aperture magnitude from the difference found between the different 
magnitude types. The aperture corrected magnitudes are presented in Table O 
Because our photometric cahbration was based on matching SExtractor's automatic aperture 
magnitudes to the z' filter photometric catalog from the COSMOS porject, this procedure 
should provide an unbiased estimate of the total (aperture-corrected) AB narrowband mag- 
nitudes for our objects. This approach also provides some robustness to crowding, thanks 
to the relatively small (1") apertures used for color measurements described in Section [3?T1 
We show our final sample in the A^i?9680 — Ks vs z' — A^i?9680 diagram represented in 
Figure [3] with the different possible contaminants described above. The colors for L- and 



T- Dwarfs have been computed using the L- an d T- Dwarfs library from lDahn et al.l ( 120021 ). 
We used GAL AXEV iBruzual fc CharlotI ( 120031 ) to model the color tracks of early and dusty 
galaxies using the Padova 1994 evolutionary tracks with a Salpeter IMF. 

We also report in Table [3] the lower lir nits of the rest-frame equiv alent widths [EW^ 
derived from the photometric data following iMalhotra fc Rhoadd (120021 ). defined as: 



^ Wrest (A) 



Ar_B9680 



AAz - /zAA 



iVS9680 



fz — fi 



iVB9680 



X 



(2) 



where fNB9680 is the observed flux in the narrow-band combined image, fz is the ob- 
served flux in the z' broad-band image, AX^smso and AXz are the width of the A^i39680 
filter (90A) and the z' band filter (928A) respectively. 

Our first object, LAE#1 in Table [31 present a detection in Ks band (see Figure S]). This 
could be due to a very red continuum slope. However, looking at the Spitzer/IRAC data in 
3.5/im and 4.8yum, none of the LAE candidates are detected. Alternatively, it could be due 
to the presence of another line, such as Mgll. Looking at the UKIRT/WFCAM2^ data in J 
band (J5o% = 25.6, AB, 5a) for all our high redshift candidates, LAEt^I is detected in this 
data with a magnitude of 24.8 (AB) and a SNR(J)~ 10. From the HST/NICMOS*^ data 
available in H band (i^5o% = 26.7, AB, 5a), LAE#3 is detected with a magnitude of 26.5 
and a SNR(H)~4. 

These detections in different broad bands confirm the reliability on these candidates. Al- 
though the other candidates have a strong single band detection, the tests realized in Para- 
graph 13. 2. II confirm that they are emission line objects. The remaining candidates are based 
on a > 5cr significant single-band detection in the narrow band filter. For Gaussian statistics, 
the false positive probability at > 5(7 is 3 x 10"'', while our survey area contains ~ 0.5 x 10^ 
independent resolution elements (based on 0.65" seeing). The number of noise spikes en- 
tering the sample should thus be ~ 1.5, comparable to the expected number of foreground 
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Fig. 3. — Color-Color Diagram z'-N B9680 vs A^i?9680-Ks showing our criteria magNBmm ~ 
magz' < —0.75mag (short dashed line) and magNBmm ~ f^dgKs < Omag (long dashed line). 
The empty triangle s ymbol correspond t o the LAE candidates. The empty squares represent 
the L-T-Dwarfs from lDahn et al.l (|2002| ). Th e green track represents ea rly-type galaxies and 
the blue track dusty galaxies generated with iBruzual fc CharlotI (120031 ) . 
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emitters. Non-Gaussian noise could increase this number, but the absence of detections in 
the negative-image test (see Section I3.2.ip support the conclusion that Sa noise spikes are 
not a major contaminant of our sample. 

Three out of six of our LAE photometric candidates are not detected in the z' band and 
we therefore use the detection limit in this band, deriving in turn lower EW limits. 



Discussion 



Variance. Two sources of variance can be involved in a such high redshift study : the 
Poisson variance and the fluctuations in the large scale distri bution of the galaxies, a lso 
called the cosmic variance. We used the on-line calculatoi0 by Trenti fc Stiavelli ( 2008 ) to 
estimate the cosmic variance. This calculator requires numerous parameters such as the area 
of the survey, the mean redshift, the reshift interval, but also the completeness value and 
several cosmological parameters. We obtained, for our sample, a value of 54% for the cosmic 
variance. The Poisson noise is estimated to 58%. The 54% uncertainty from the cosmic 
variance and the 58% from Poisson statistics are comparable and we therefore consider both 
in our error estimation. 

Luminosity Function. Our A^i?9680 filter being quite narrow (FWHM~90A) , we assume 
that the narrow-band flux is entirely coming from the Lya line. We fit to the Lya luminosity 
function of this z~6.96 LAE sample, a Schechter function, ^{L), given by 



$(L)dL = ^*{^ 



L\ dL 



' L* L* 



(3) 



in order to compare with previous high redshift works (lOuchi et al. 



Tilvi et al.ll2010l : lOuchi et al.ll2008l : lOta et al.ll2008l : iKashikawa et al. 



2010 



2006 



Hibon et al. 


2010 




Malhotra & Rhoads 



20041 ). Considering the small number of candidates in our sample, we choose to fit two out 
of three of the Schechter function parameters. We set the faint end slope of the luminosity, 
a, to a = —1.5, and determine $* and L* by minimization. We decide to obtain a best-fit 
Schechter functions for the z~7 cumulative luminosity function, which has been derived by 
considering only our photometric candidates. 

In Figure O we present the Lya Luminosity Function not corrected for detection incomplete- 
ness as a black solid line. We then use the completeness result of Section 12.2.11 to correct 
the luminosity of our objects, and find a new Schechter fit, presented as a red solid line in 
Figure This is the z 7 Lya Luminosity Function corrected for incompleteness. 



^http: / /casa. Colorado. edu/^trenti/CosmicVariance. html 



Tabic 3: Tabic of tlic z-^ 6.96 LAE candidates. 



Name 


Ra 


Dec 


NB9680 


Error 


SNR {NB9680) 


z' 


Error 


SNR{z') 


EW- (A) 


J 


Error 


LAE#1 


10:00:46.846 


02:10:16.01 


24.1 


0.2 


7.3 


26.0 


0.16 


6.6 


~41 


24.8*^ 


0.15 


LAE#2 


10:00:39.104 


02:03:02.55 


24.1 


0.18 


7.2 


>26.4 






>61 


>24.8 




LAE#3 


09:59:50.991 


02:12:19.13 


24.2 


0.2 


6.4 


>26.4 




2 


>47 


>24.8 


... 1 


LAE#4 


10:00:58.529 


02:12:56.51 


24.2 


0.18 


6.1 


>26.4 






>47 


>24.8 


... ^ 


LAE#5 


10:00:42.528 


02:11:30.39 


24.3 


0.19 


6.1 


>26.4 




0.2 


>48 


>24.8 


1 


LAE#6 


10:00:37.940 


02:11:59.28 


24.4 


0.2 


5.6 


>26.4 






>41 


>24.8 





^ In the restframe 



^ from UKIRT/WFCAM2 data 
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Fig. 4. — Thumbnail images of all candidates listed in Table |3l Each window is 15" x 
15". Objects names and passbands are located above and to the left of the thumbnails, 
respectively. 
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Fig. 5. — Best-fit Schechter function for the cumulative z~7 Lya luminosity function. The 
red solid line is the best-fit z~7 LF with completeness correction. Our candidates are rep- 
resented as triangles (plain: not corrected for completeness, empty : with completeness 
correction). The error bars represen t the Poisson er rors and are i dentical fo r plain and 
empty points. The z~6.96 LAE from llye et al.l (120061 ) is the square, lOta et al.l (12010) can- 
didates are pentagons. Also represented here is the z=6.5 Lya LF from lOuchi et al.l (120101 ) 
and the z=5.7 Lya LF from lOuchi et al.l (120081 ). 
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Sample Incompleteness. We create a grid pattern of 15000 objects on a mock image, an- 
drun SExtractor for different pfiotometric apertures in double image mode, using the ^r'-band 
image as the measurement image. We remark that by increasing the photometric aperture 
size, the number of objects matching the optical criterion (Criterion ^2 in Table [2]) decreases. 
For an aperture of 5 pixels, we recover 75%±0.65%, for 10 pixels 63%±0.6% and for 20 pixels 
46%±0.5%. We choose a photometric aperture of 5 pixels for the objects catalogs we used 
for the high redshift LAEs selection. We know then that we could miss 25% of the objects 
due to the photometric aperture size we choose. This corresponds to the possibihty that we 
missed ~1.5 objects in our high redshift sample. Assuming one more object in our sample, 
our conclusion about the best-fit LP will not change. 



Interpretation. The previous studies presenting z~7 Lya emitters (lOta et al.l 12008 



20101) has lead to the first spectroscopic confirmed z~7 LAE, called lOKl. Their survey cov- 
ers an area of 876 square arcmin with the filter NB973 (AA = 200A, Ac = 9755A) and reaches 
a 50% completeness of NB973=25.6 (AB, 5a") (equivalent to a flux of 1.36e~^^ergs~^, cm~^) 
in the SXDS (Subaru/XMM-Newton Deep Survey) and NB973=25.3 (AB, 5a) (equivalent 
to a flux of 1.8e~^^ergs~^, cm~^) in the S DF (Subar u Dee p Field). IOK-1 has a flux of 



2 X 10~^^ergs~^ cm~^. From the Table 1 of lOta et al.l (120081 ) . we are able to obtain a lower 
limit for the rest-frame equivalent width of lOKl, EWrest ~ 49A, using Equation [2l This 
rest-frame equivalent width is in agreement with the rest-frame equivalent l imit l ower limit 



we found for our candidate sample and presented in Table Ei lOta et al.l (120101 ) find four 
new photometric candidates in a survey covering the Subaru/XMM Newton Deep Survey 
Field with Suprime-Cam and reaching a depth limit of NB973=25.4, corresponding to 72% 
completeness. These candidates are represented by pentagons in the Figure [51 



From lYoshida et al.l (120061 ) . lOta et al.l (120081 ) estimated therefore a possible z~7 Lya 



Table 4: Best fit Schechter LF parameters for a = —1.5 



Redshift log(L*(ergs~^)) log(<l>*(Mpc~^)) 



6.96^ 


42.56+IJi 


-2.0lt'of 


6.96(2) 


42.8l°:J^ 


-3.441°:?° 


6.5(3) 


42.64+°-J 


-3.07l°:i^ 


5.7(4) 


42.8t',il 


q 1 1 +0.29 
"-'•-'--'--0.31 



Refer ences. (1) d e rived for our sample; (2 lOta et al.l ( 120081 ): (3) 
Ouchi et all toO^ 



Ouchi et al, 



(120101); (4) 
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LF with a pure lumino sity evolution of L*^j = 0.58L*^5 7, with -f/^=5.7 = 1.08.10^^ergs~^ 
(jShimasaku et al.ll2006l ). This inferred z~7 Lya LF predicts fewer LAEs than seen in our 
photometric candidate sample. Confirni ation of 1-2 of our candidates would approximately 
match the prediction in lOta et al.l (|2008[ ) and would modestly exceed their measured number 
density. 



We show in Figure the cumulative z~7 LAEs LF obtained after correcting our points 
for the aperture and the detection completeness. This completeness correction has been 
applied by number weighting according to the NB9860 magnitude. The best-fit parameters 
do not vary significantly before and after correcting from the completeness, as seen in Figure [5] 
between the black solid line and the red solid line (without and with the completeness 
correction, respectively). 

By considering only our photometric candidate sample, we do not observe any strong L* 
or $* evolution between z=5.7 and z~7, and therefore contradict a possible L* evolution 
between z=6.5 and z~7. If none of our candidates is a real z~7 LAE, we can then put an 
upper limit on the z~7 Lya LF, which will help constrain the neutral fraction of the IGM. 



5. Conclusions 

We observed 465 arcmin^ from the COSMOS field using the narrow-band imaging tech- 
nique on Magellan/IMACS with the A^i?9680 filter, in order to target the Lya line at z~6.96. 
We obtained a comoving volume of ~ 72000Mj9c^. After applying our selection criteria and 
verifying that our selection was not contaminated by low-redshift emitters and transient 
objects, we obtain a sample of six z~6.96 LAEs. From this photometric sample, we are able 
to infer a possible z~6.96 Lya Luminosity Function. We find no evolution in luminosity 
function from z=6.5 to z~6.96, if a majority of our sources are confirmed. 
It is now crucial to obtain spectroscopic follow-up observations to reveal the real nature of 
these objects and establish a firm conclusion on the z~6.96 Lya Luminosity Function. 



The authors would like to thank the referee for very useful comments. Las Campanas 
Observatory staff, the IMACS team and the National Science Foundation. 
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